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TIME DELAYS IN THE GRAVITATIONALLY LENSED QUASAR H1413+117 (CLOVERLEAF) 
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ABSTRACT 

The quadruple quasar H1413+117 (z s = 2.56) has been monitored with the 2.0 m Liverpool Tele- 
scope in the r Sloan band from 2008 February to July. This optical follow-up leads to accurate light 
curves of the four quasar images (A-D), which are defined by 33 epochs of observation and an average 
photometric error of ~ 15 mmag. We then use the observed (intrinsic) variations of ~ 50—100 mmag 
to measure the three time delays for the lens system for the first time (la confidence intervals): Atab 
= —17 ± 3, Atac = — 20 ± 4, and Atad = 23 ± 4 days (Ar^ = Tj — r^; B and C are leading, while 
D is trailing). Although time delays for lens systems are often used to obtain the Hubble constant 
(Ho), the unavailability of the spectroscopic lens redshift (zi) in the system H1413+117 prevents a 
determination of Hq from the measured delays. In this paper, the new time delay constraints and 
a concordance expansion rate (Hq = 70 km s _1 Mpc -1 ) allow us to improve the lens model and to 
estimate the previously unknown zi. Our la estimate zi — 1.88^q is an example of how to infer the 
redshift of very distant galaxies via gravitational lensing. 
Subject headings: gravitational lensing — quasars: individual: H1413+117 



1. INTRODUCTION 

The time delay between two images of a gravitation- 
ally lensed source depends on the distribution of mass 
in the lens and the c u rrent expansion rate of the Uni- 
verse (jRefsdal I fl96i [1966). This expansion rate is 
quantified by the Hubble constant Hq. If the source 
is a quasar, the intrinsic quasar variability may be 
used to determine the time delays between its multi- 
ple images. Thus, observed delays for lensed quasars 
lead to valuable information on Hq, provided lensing 
mass dist ributions can be constrained by observational 
data fe.g..lKochanek fc Schechterl|2004l:lSchechter 112004 
iSaha et al. 1 120061: 1 Jackson Il2007t lOguri 1 120071) . Future 
large samples of lens systems could be useful tools to 
obtain accurat e estimates of the m ain cosmological pa- 
rameters (e.g.. lDobke et al."ll2009f ). 

For a given lensed quasar, each time delay between 
two of its images is indeed scaled by a factor containing 
Hq, the lens redshift z\, and additional physical param- 
eters. If z\ is known (this is the usual situation), the 
Hubble constant and the lensing mass distribution can 
be simultaneously deduced by using a parametric lens 
scenario and a set of observational cons traints (including 
inform ation on the time delay(s); e.g.. iSchneider et al. I 
2006). However, a large set of constraints is required to 
reliably determine both cosmological and galactic prop- 
erties. 

Accurate observations of a quadruply imaged quasar 
in a well-studied galaxy field bring an excellent oppor- 
tunity to study in detail Hq and the mass distribution of 
the gravitational lens. Apart from the three independent 
time delays, and the positions and fluxes of the four im- 
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ages, data on the neighbour galax ies also constraint the 
lens scenario (e.g.. Uack son 2007). Alternatively, if zi is 
very uncertain or unknown, one may infer a lens redshift 
value from the time delay measurements and the rest of 
constraints (using a value of Hq from other experiments). 

In this paper we f ocus on the quadru ple quasar 
H1413+117 fCloverleaf: iMagain et al. Ill988l). lyin g at a 
redshift z s = 2.558 (e.g.. iBarbainis et al. Tll997T l. Al- 
though the lens redshift is curren tly unknown, severa l 
neighbour galaxies w ere detected bv lKneib et all ([1998). 
iKneib et al. I (1998) found the main lensing galaxy Gl 
amid the four quasar images, and presented astrometric 
and photometric data for additional galaxies surrounding 
H1413+117. They analysed a secondary lensing galaxy 
(G2) close to Gl, as well as some other objects proba- 
bly belonging to an overdensity (galaxy group / cluste r) at 
photometric redshift z ove ~ 0.9. iFaure et al. I (|2004D also 
found evidence for the presence of two different overden- 
sities at z ove = 0.8 ± 0.3 (corresp onding to the structure 
discovered bv lKneib et al. Ifl998f ) and z ove — 1.75 ± 0.2, 
which could contribute noticeably to the lensing poten- 
tial. 

Very recently, iMacLeod et al. I (|2009D have used the 
positions of Gl, G2, and other candidate lensing galax - 
ies, the quasar image positions (jTurnshek et al. I I 1997). 
new mid-IR flux ratios, and some priors to constrain the 
lensing mass distribution. They have concluded that the 
galaxy pair G1-G2 and an external shear (likely related 
to the observed galaxy overdensities) are required to ex- 
plain the observati ons. From an opt i cal m onitoring over 
the 1987-94 period. lOstensen et al. I (|1997j ) also reported 
a quasi-simultaneous variability of the four quasar im- 
ages. Unfortunately, the scarce sampling did not permit 
them to estimate the time delays in the system. Hence, 
the measurement of time delays for the Cloverleaf quasar 
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should improve knowledge about the lens (mass and red- 
shift), as well as allow completion of new studies (e.g., 
microlensing variability) . 

In Section 2 we present new optical light curves of 
H1413+117. In Section 3, from these light curves, we 
estimate the time delays between quasar images. In Sec- 
tion 4 we compare the most recent lens scenario with all 
relevant observations. In Section 5 we discuss our results 
and put them into perspective. 

2. OBSERVATIONS AND LIGHT CURVES 

We observed H1413+117 from early February to late 
July of 2008, i.e., for 6 months. All observations were 
made with the 2.0 m fully robotic Liverpool Telescope 
(LT) at the Roque de los Muchachos Observatory, Ca- 
nary Islands (Spain), using the RATCam optical CCD 
camera (binning 2x2). The global database consists of 
61 exposures of 300 s in the r Sloan filter 1 . These original 
exposures (frames) are almost regularly distributed over 
the whole observation period (there is only a significant 
21-day gap in 2008 February), with an average sampling 
rate of about one frame each three days. In Fig. 1, a 
combined LT frame shows the blended quasar images (at 
the centre of the field), two relevant stars (the control 
star S40 and the PSF star S45 , which correspond t o the 
objects 40 and 45 in Fig. 1 of iKavser et al. IH990f ). and 
other relatively bright objects. The left bottom corner 
of Fig. 1 illustrates the positions (crosses) and names 
(A-D) of the four quasar images. 

The four quasar images are separated by ~ 1" (see Fig. 
1), so we only consider 33 high-quality frames to make 
the light curves (photometry on the other 28 frames is 
discussed in the last paragraph of Section 3). This selec- 
tion is based on the FWHM of the seeing disc measured 
in each frame (FWHM < 1'.'5), as well as the signal-to- 
noise ratio (SNR) of the 18.16-mag control star (SNR 
> 150). We note that this star (S40 in Fig. 1) has an 
r-band brightness similar to those of the quasar images 
(~ 17.9-18.4 mag). Each SNR value is calculated within 
an aperture with radius equal to the frame FWHM. 

We determine the instrumental fluxes of the four 
quasar images through point-spread function (PSF) fit- 
ting. As the main l ensing galaxy is ver y faint in the r 
band (R > 22.7 mag jKneib et al. |[l998f) . our photomet- 
ric model includes four stellar-like sources (i.e., 4 empir- 
ical PSFs) plus a constant background. The empirical 
PSF is derived from the 16.69-mag star in the vicinity 
of the lens system (S45 in Fig. 1). In order to obtain 
accurate and reliab le fluxes, we use the well-tested IM- 
FITFITS software (TMcLeod et al. II1998D . incorporating 
th e relative positions of t he B-D images (with respect to 
A: lTurnshek et alHll997h as constraints. Thus, the code 
is applied to all selected frames (see above), allowing 7 
parameters to be free. These free parameters are the 
position of the A image, the fluxes of the A-D images, 
and the sky brightness. We also infer SDSS magnitudes 
from the relative instrumental magnitudes of the lensed 
quasar. The S45 star is taken as reference for differential 
photometry. In Table [T] we present the r-SDSS magni- 
tudes of the four images and the control (S40) star. 

1 The pre— processed frames are publicly available 
on the Liverpool Quasar Lens Monitoring archive at 
http : //dc .zah.uni--heidelberg. de/liverpool/res/rawf rames/q/f 



From the standard deviation of the magnitudes of 
the S40 star over the whole monitoring period, we ob- 
tain a typical error of 0.006 mag. This global scatter 
agrees with the standard deviation between magnitudes 
on consecutive nights divided by the square root of 2, 
as expected on theoretical grounds. To estimate typical 
photometric errors in the quasar light curves, we then 
use the standard deviations between magnitudes having 
time separations < 1.5 days (true variability is negligible 
on this very short timescale), which are divided by the 
square root of 2. The resulting uncertainties are 0.010 
(A), 0.012 (B), 0.018 (C) and 0.018 (D) mag. As a sum- 
mary, we achieve ~ 1-2% photometry and reasonable 
sampling rate (~ 6 data per month). Moreover, the light 
curves of the A-D images show significant variations of ~ 

0. 05-0.1 mag (see Fig. 2). For example, the almost par- 
allel fading by ~ 0.1 mag of A-D (over the last 100 days 
in Fig. 2) suggests intrinsic variability. This is promising 
to derive time delays. 

3. TIME DELAYS 

There are four different image ray paths for the Clover- 
leaf quasar, so travelt i me (r ) varies from image to im- 
age ([Schneide r et al. I Il992f ). Assuming that the ob- 
served magnitude fluctuations are basically originated 
in the source quasar (intrinsic variability), we use two 
well-known cross-correlation techniques to measure time 
delays between quasar images Ar,j = Tj — Ti, where 

1. j = A, B,C, D. These techniques are the dispersion 
(D 2 ) and reduc e d chi-square (\ 2 ) minimizations (e.g., 
iPelt et alJ [19961 lUllan et al. 1 12006T ). Despite the ex- 
isten ce of other methods for determining time delays 
Ce-g- lKundic et al. Ifl997l : I Gil-Merino et al. 1 120021. most 
methods work in a similar way, and in most cases one 
does not need to carry out an exhaustive analysis. After 
deriving delays we discuss the intrinsic variability hy- 
pothesis at the end of this section. 

We focus on the AB, AC, and AD comparisons, i.e., the 
A light curve is compared to the other three brightness 
records (B-D). The D 2 and x 2 minimizations are char- 
acterized by a decorrelation length (<5) and a bin semi- 
size (a), respectively. To simultaneously avoid very noisy 
trends and loss of signal (due to excessive smoothing), we 
take a = S = 15 days. For a = S = 15 days, the spec- 
tra between —75 and +75 days include global and local 
minima. However, these local minima do not play an im- 
portant role in the estimation of the time delays Atab , 
Atac, and At ad (see below). 

For a given cross-correlation method, we follow two 
different approaches to generate synthetic light curves 
and determine time delay errors. In the first approach 
(which is called NORMAL), we do not make any hypoth- 
esis on the underlying signal, but the observational noises 
(in the four light curves) are assumed to be normally dis- 
tributed. Therefore, one obtains a synthetic light curve 
of an image by adding a random quantity to each bright- 
ness in the observed record. These random quantities are 
realizations of a normal distribution around zero, with a 
standard deviation equal to the standard deviation be- 
tween observed magnitudes on consecutive nights. We 
produce 1000 synthetic light curves of each image, and 
thus, obtain 1000 delay values for each pair (AB, AC, 
and AD) and the corresponding 68% confidence inter- 
ormvals. In the second approach, we use a bootstrap pro- 



Time Delays in H1413+117 



3 



cedure (BOOTSTRAP; e.g.. lEfron fc Tibshirani I [l993h . 
First, we derive a combined light curve from the best so- 
lutions of Atab, At^c, and At ad (global minima of the 
three spectra), and the associated magnitude differences. 
This combined curve is then smoothed by a 20-day filter. 
Second, the combined and smoothed curve is assumed to 
be a rough reconstruction of the underlying signal, and 
thus, the residuals for each image are taken as sets of 
errors. These four sets are resampled to infer 1000 boot- 
strap simulations (synthetic light curves) of each image. 
To measure the time delays (68% confidence intervals), 
we compute 1000 delay values for each pair of images and 
analyse the three distributions. 

The time delay measurements are presented in Table [2j 
For each pair of images, the four measurements are con- 
sistent with each other. Table [2] also indicates the exis- 
tence of an average offset of about 3 days between the 
results from the NORMAL approach and those derived 
from BOOTSTRAP simulations. This offset may be due 
to slightly biased reconstructions of the underlying sig- 
nal (BOOTSTRAP procedure). We adopt At AB = -17 
± 3, Atac = -20 ± 4, and At A d = 23 ± 4 days as 
our final la measurements (x 2 and NORMAL). The cor- 
responding r-band magnitude differences are Aitiab = 
0.155 ± 0.006, Am AC = 0.322 ± 0.011, and Am A D = 
0.501 ± 0.007 mag (la confidence intervals). These give 
the magnification (flux) ratios: B/A = 0.867 ± 0.005, 
C/A = 0.743 ± 0.007, and D/A = 0.630 ± 0.004. 

From the central values in the time delay and mag- 
nitude difference intervals, one can obtain a final com- 
bined light curve, i.e., the A light curve together with 
the magnitude- and time-shifted B-D records. In Fig. 
3 we display this combined light curve (filled symbols). 
Using a 20-day filter (filter semisize = 10 days), a pos- 
sible reconstruction of the underlying intrinsic signal is 
also drawn in Fig. 3 (solid line). The standard deviation 
between the magnitudes in the combined record and the 
reconstruction is about 0.015 mag (see the shaded area 
in Fig. 3, which represents 0.015-mag deviations from 
the reconstruction), in good agreement with the average 
photometric error in the quasar light curves (see the error 
bar in the lower left corner of Fig. 3). This suggests the 
absence of significant microlcnsing signals, and strength- 
ens our hypothesis that observed variat ions mainly have 
an int rinsic origin. We also note that lOstensen et al. I 
(|1997t ) reported on an almost parallel variation in bright- 
ness (R band) of the four quasar images over a 7-year 
period (see the middle panels in Figs. 2-3 of that pa- 
per). Moreover, our 6-month monitoring period was long 
enough to find significant variability and determine time 
delays, but it was not long enough to detect substan- 
tial microlcnsing variations. Typical microlensing gradi- 
ents of ~ 10~ 4 mag day"" 1 are expected in the rR bands 
/e.g.. iGavnullina et al. 1 120051 : iFohlmeister et al. 1 120071 : 
iShalvapin et al. 112009ft . 

Once the time delays are measured and the combined 
light curve is drawn in Fig. 3 (filled symbols), we can 
discuss the accuracy of the magnitudes derived from the 
poor-quality frames (see Section 2). These 28 exposures 
with FWHM > l'/5 and/or SNR < 150 seem to be use- 
less, since they lead to a very noisy combined record and 
confusion. One frame has a very poor image quality, 
so we do not extract quasar magnitudes. Some of the 
magnitudes associated with the other 27 poor-quality 



frames are shown in Fig. 3 (open symbols). The rest are 
extremely noisy, and their values are outside the magni- 
tude range in Fig. 3. In brief, 28 out of the 61 original 
frames have either an excessive blurring, or an insuffi- 
cient signal, or both of them, so they do not produce 
accurate quasar light curves. 

4. IMPROVED LENS MODEL AND LENS REDSHIFT 

In a model-independent way, the image and main lens 
positions for the H1413+117 system are useful to deter- 
mine the ordering of the time delays (jSaha fc Williams! 
2003). However, the quasar images B and C (associated 
with minimum arrival times) are almost equidistant from 
the main lens, so it is difficult to distinguish the leading 
image in this system. In any case, intrinsic variations 
should be firstly observed in light curves of these two im- 
ages, and later in records of A and D. While D is the trail- 
ing image (it is the closest to the main lens), A should be 
characterized by an intermediate arrival time. Our time 
delay measurements agree with this time-ordering of the 
images. Detailed lens models predict that C is leading 
(jChae fc Turnshek iflOQl IMacLeod et al. ll27)09ft . How- 
ever, we cannot confirm this prediction at la confidence 
level, since the Atab and Atac intervals overlap each 
other (a direct measure Atbc is neither useful to decide 
o n the leading image). 

IMacLeod et al. I ([2009D reported how a relatively sim- 
ple lens model reproduces the observed positions and 
mid-IR flux ratios of the four quasar images. These mid- 
IR flux ratios are insensitive to extinction (long wave- 
length) and microlcnsing (large emission region). The 
MacLeod et al.'s main solution (see the second column 
in Table [3]) relies on an observationally motivated sce- 
nario. This consists of a background point-like source 
(quasar) that is lensed by a singular isothermal ellipsoid 
(main lensing galaxy Gl), a singular isothermal sphere 
(secondary lensing galaxy G2), and an e xternal shear 
(likely produced by gala xy overdensities; iKneib et al. I 
119981: iFaure et al.ll2004D . Although the position of Gl 
was constrained by observations, it was allowed to vary 
during the fitting procedure. The singular isothermal 
sphere was placed at the observed position of G2, and 
MacLeod et al. also assumed priors on the ellipticity of 
Gl (eci = 0.0 ± 0.5) and the strength of the external 
shear ( lext = 0.05 ± 0.05). 

Here, we use the MacLeod et al.'s lens scenario and a 
concordance cosmolog y: H n = 70 km s" 1 Mpc" 1 , fl m 
= 0.3, fl A = 0.7 (e.g.. iSpergel et al. Il2003l ). The theo- 
retical time delay between two lensed images includes a 
cosmological scale factor T = (1 + zi)(DiD s /cDi s ), where 
c is the velocit y of light and D denote s angular diameter 
distance (e.g.. ISchneider et al. 111992ft . The angular di- 
ameter distances are determined by the cosmology, and 
the lens and source redshifts, z\ and z s . Thus, using a 
concordance cosmology and the observed reds hift of the 
source z s — 2.558 fe.g. JBarbainis et al. 111997ft . the scale 
factor exclusively depends on z\. Our goal is to find a 
good fit to all observations of interest, i.e., the image po- 
sitions and (mid-IR) fluxes, and the three time delays 
in Section 3. We also use the constrai nts and priors on 
the G 1-G2 positions, eat, and ^ ext by IMacLeod et al71 
(2009). With respect to the MacLeod et al.'s framework, 
we add 3 new observational constraints (time delays) and 
one new free parameter (z;), so the degrees of freedom 
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(dof) change from 5 to 7. 

Th rough the GRAVLENS 2 software package (jKeeton I 
120011 ) , we find a solution with x 2 /dof = 7.5/7. Our main 
results are shown in the third column of Table [3] Most 
lensing mass parameters of this acceptable solution are 
close (deviations less than 10%) to those in the second 
column of the same table (see al so the third column in 
Table 3 of lMacLeod et al. 112009th However, it is inter- 
esting to note that the new external shear strength (j ex t 
= 0.11) is slightly larger than the previous one, whereas 
the mass scale of G2 (bo2 — 0('45) is slightly smaller (de- 
viations of 25-30%). On the other hand, the best-fit lens 
redshift indicates the existence of a very distant galaxy 
pair G1-G2. In Fig. 4 we draw the x 2 — zi relation- 
ship, which permits us to obtain confidence intervals for 
zi. The ler determination is zi = 1.88^° (dark shaded 
area in Fig. 4). Fig. 4 also shows the 2a (95%) confi- 
dence interval: 1.65 < zi < 2.03 (whole shaded area). 

5. DISCUSSION 

In the lens system H1413+117, the main l ensing galaxy 
Gl is a very faint object (i? > 22.7 mag; IKneib et al. I 
1998), surrounded by four close and relatively bright 
quasar images (R ~ 18 mag). Thus, it is very difficult 
to separate the spectrum of the galaxy from those of the 
quasar images. Moreover, the quasar spectra indicate the 
presence of intervening objects (absorptio n lines) at dif- 
feren t redshifts less than z s = 2.56 fe.g.. IMonier et aT~l 
119981 and references therein), so one cannot decide on 
the redshift of Gl. The secondary lensing galaxy G2 and 
other galaxies in the vicinity of t he lensed quasar ar e also 
very faint objects (R > 23 mag; Kn eib et al. Ill998f ). and 
spectroscopic redshifts of these galaxies are not available 
yet. Photometric data of field galaxies are consistent 
with the presence of two galaxy ove rdensities at z ?ve = 
0.8 ± 0.3 and z m )e . = 1.75 ± 0.2 (IKneib et al. I [19911 
iFaure et al. I [20041 . For example, the galaxy G2 ha s 
a photometric redshift of about 2 (jKneib et al. Ill998f ). 
Our gravitational lensing estimate of the redshift of Gl- 
G2: zi = 1.88^0 11 (l* 7 interval), is in reasonable agree- 
ment with the photometric redshift of G2 and the most 
distant overdensity, as well as the absorption system at 
z a bs = 1-87. However, the nearest group/cluster is far 
away from the principal gravitational deflector (galaxy 
pairGl-G2). 

So far it is not considered a possible uniform external 
convergence, i.e., n ext — 0. If n ext ^ 0, then fits to K ext = 
lens scenarios (lens models with n ext — 0) must be con- 
veniently rescaled. Original estimates of Ho (if that were 
the case) or z\ (in our case) also require suitable correc- 
tions. This is because of the so-called mass sheet degen- 



Falco et al. 1 119851: iGqrenstein et al~1 119881: 
Nakaiima et al. Il2009f) . For H1413+117. it 



eracy (e.g., 
[Sal^[2000l : 

is unclear what is the main perturber producing external 
shear, and perhaps external convergence. For example, 
infrared photometry of the neighbou r galaxy H2 is con - 
sistent with a redshift of about 2 (jKneib et al. i ll 998). 
so it could be located in the principal lens plane. More- 
over, the orientation of the external shear, (9 7cxt = 45?4, 
is in the direction of this galaxy (e.g., see Fig. 4 of 
iMacLeod et~al~l 12009). Thus, H2 and some other re- 

2 http: //redf ive .rutgers . edu/$\sim$keeton/gravlens 



lated objects (belonging to the very distant overdensity) 
might produce most of the external shear and a negligi- 
ble convergence. Apart from this optimistic perspective, 
one may also consider that the main perturber is the 
most distant group/cluster as a whole. If this were true, 
the external convergence at the quasar positions would 
be K ex t = j ex t ~ 0.1 for a singular isothermal sphere. 
Taking n ext ~ 0.1 into account, the derived lens redshift 
increases by ~ 3% (+ 0.05). This increase represents 
only one half of our \a uncertainty in zi (~ 0.10; see 
above) . 

We can also quantify the gravitational influence of the 
nearest overdensity. A weak lensing analysis of this struc- 
ture indicated a sh ear direction of 45° (see the last col- 
umn m Table 5 of IFaure et aFI I2004D . This shear di- 
rection coincides with our strong lensing determination 
in Table [H which suggests that the nearest overden- 
sity might pl a y a no ticeble role in the lens phenomenon. 
IFaure et al. I (|2004l ) also determined an u pper limit on 
the s h ear strength: 7 < 0. 17 (see Fig. 7 of IFaure et al. I 
[2001 . iKeeton 1(120031 ) and lMomcheva et al. I (|2006l) dis- 



cussed the effective convergence and shear when a per- 



turber does not lie at 



zu e.g. 



^per 



< zi . Assuming 



a singular isothermal sphere to describe the perturber, 



ieff = Jeff ~ (1 - 13) J, where (3 ~ if z v 



zi and 



> if z per < zi (see Appendix of iMomcheva et a"l~l 
2006). From the involved redshifts and the upper limit 
on the shear strength, we infer that K e ff = Jeff < 0.03. 
Thus, the nearest structure cannot account for the ex- 
ternal shear in the lens system, and it likely generates 
no more than 10% of the total shear strength. The con- 
straint on the effective convergence leads to a negligible 
increase in z;, i.e., Az; is below + 0.01. In some other lens 
systems, there are also perturbers at redshifts less than 
those of the principal lensing objects, which produce a 
few hundredths of external convergence and shea r (e.g., 
iFassnacht fc Lubin 112001 iFassnacht et al. II2006D . 

Our lens scenario incorporates singular isothermal 
mass profiles. Despite the fact that these profiles lead 
to an acceptable fit (x 2 /dof = 7.5/7), other distribu- 
tions of mass (including a core, some deviation from the 
isothermal behaviour or both ingredients) could also lead 
to good fits for the observation al data. This is the well- 
known profile degeneracy (e.g., Uackson |[2007l ). An ex- 
haustive study of mass distributions consistent with the 
available observational constraints is out of the scope of 
this paper. However, we check the influence of non- 
isothermal profiles of Gl on z\ estimations. The power- 
law index of Gl (a) is assumed to be around 1 (isother- 
mal index), and additional fits with singular a/1 pro- 
files are done. For a = 1.1, the solution is characterized 
by Ax 2 = X 2 -X 2 (« = 1) = 0.1. For a = 0.8-0.9, we find 
a very modest improvement in x 2 , Ax 2 = —0.1. The full 
range a = 0.8-1.1 leads to best -fit values of zi within 
our la "isothermal" estimate in Section 4. 

Microlensing variations in H1413+117 may shed light 
on the nature and s t ructure of the source quasar (e.g., 
iLewis fc Bellel [19981 iPopovic et al. I [20061 ). The three 
time delay measurements in Section 3 and the new data 
on the lens (mass and redshift) in Section 4 are use- 
ful tools for analyses of microlensing variability. Once 
the delays are known, it is possible to properly com- 
pare quasar light curves and to search for microlensing 
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signals. Moreover, the (lens and source) redshifts and 
the improved lens model allow construction of microlens- 
ing magnificatio n patterns and simu lated microlensing 
light curves (e.g.. iWambsganss Hl990T ). The key parame- 
ters for microlensing simulations are the convergence and 
shear strength at the positions of the quasar images. To 
address the space distribution of both convergence and 
shear, we consider our lens model in the third column 
of Tabled This gives (k a , -/ a ) = (0.51, 0.58), (k b , 
7b) = (0.52, 0.32), (k c , 7c) = (0.48, 0.37), and (k d , 
Id) = (0.57, 0.65). Although we cannot rule out the 
existence of an external convergence at a level of 0.1, 
this scenario is uncertain (see above). Hence, we do not 
take into account any external convergence due to galaxy 
overdensities along the line of sight to the lensed quasar. 

Our lens model causes different magnifications at the 
four positions of the quasar images. The model magnifi- 
cation ratios are B/A = 0.82, C/A = 0.78, and D/A = 
0.45. All these ratios agree, as expect e d, wi th the mid- 
IR measurements by iMacLeod et al. I ()2009f ) . However, 
the model ratios are not included in the error bars of 
our optical (r-band) flux ratios in Section 3. These dif- 
ferences between optical and model ratios are probably 
due to dust extinction and microlensing magnification. If 
microlensing is currently playing a role (e.g., the optical 



continuum of the image D cou ld be magnified by mi- 
crolensing; lAnguita et al. Il2008l . and references therein) , 
it should be a long-term effect that induces small (opti- 
cal) flux variations on time scales of several months. This 
microlensing variability scenario is supported by recent 
studi es for other systems with non-local lens galaxies 
(e.g.. iGavnullina et al. 1 120051: iFohlmeister et al. I [20071 
IShalvapin et al. 1120091 ) 
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We thank the Liverpool Telescope staff for kind interac- 
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We use datasets taken from the Sloan Digital Sky Survey 
(SDSS) Web site, and we are grateful to the SDSS team 
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supported by the Spanish Department of Education and 
Science grant AYA2007-67342-C03-02, and University of 
Cantabria funds. 
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FlG. 1. — r— band Liverpool Telescope imaging of the Cloverleaf quasar. We combine the 5 best frames in terms of seeing (total exposure 
time = 1500 s, FWHM = 0'.'87). We then zoom into the central region and perform linear interpolation between adjacent pixels (box in 
the left bottom corner). The four quasar images and two relevant stars are properly labeled. 
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Fig. 2. — Light curves of four quasar images A-D and the control star S40. The stellar record is shifted by +0.3 mag to facilitate 
comparison. 
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Fig. 3. — Final combined light curve of H1413+117. This is made from the A light curve and the magnitude- and time-shifted B-D 
records (filled symbols). To create the combined record, we use the time delays Atab = — 17, Atac = — 20, and Atad = 23 days 
(shifts in time), and the magnitude differences Atuab = 0.155, Arrive = 0.322, and Atuad = 0.501 mag (shifts in magnitude). We also 
display a possible reconstruction of the quasar signal (solid line), as well as the standard deviation between the combined record and the 
reconstruction (shaded area). The error bar in the lower left corner represents the average photometric error in the quasar light curves. To 
check the accuracy of the quasar magnitudes in 27 poor-quality frames that are not used in the time delay analysis, we also compare some 
of these additional magnitudes (open symbols) with those derived from the high— quality exposures (filled symbols). Several additional data 
(for the poor— quality frames) are outside the magnitude range from 18.0 to 17.8. 
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Lens redshift 



Fig. 4. — Estimation of the previously unknown lens redshift. We show both the lcr (dark shaded area) and 2cr (whole shaded area) 
confidence intervals. We estimate z\ via gravitational lensing, using observational constraints on the lensed quasar images (positions, fluxes, 
and time delays), as well as data on the neighbour galaxies. 
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TABLE 1 
Photometry of H1413+117. 



Til 1 TT1 1 IPf 


rivil Hfitp a 


MJD-54000 


FWHM b 


SNR C 


A d 


B° 


C f 


D g 


S40 h 


1 


Feb 1 


498.1885 


1.08 


235 


17.884 


18.077 


18.214 


18.384 


18.168 


2 


Feb 4 


501.1370 


1.07 


232 


17.886 


18.084 


18.205 


18.375 


18.162 


3 


Feb 5 


502.1893 


1.26 


233 


17.867 


18.083 


18.235 


18.405 


18.167 


4 


Feb 29 


526.0175 


1.45 


164 


17.935 


18.052 


18.149 


18.384 


18.163 


5 


Mar 1 


527.0400 


1.24 


186 


17.932 


18.033 


18.162 


18.369 


18.160 


6 


Mar 3 


529.0227 


1.44 


173 


17.901 


18.029 


18.171 


18.396 


18.160 


7 


Mar 7 


533.0942 


1.19 


189 


17.906 


18.041 


18.150 


18.409 


18.171 


8 


Apr 12 


569.1711 


1.31 


192 


17.872 


18.028 


18.213 


18.337 


18.165 


9 


Apr 22 


578.9261 


1.02 


192 


17.870 


18.015 


18.166 


18.353 


18.164 


10 


Apr 24 


580.9120 


1.32 


202 


17.864 


18.002 


18.187 


18.364 


18.163 


11 


May 1 


587.9478 


1.03 


212 


17.847 


18.003 


18.190 


18.351 


18.161 


12 


May 3 


589.9377 


1.22 


210 


17.877 


17.989 


18.173 


18.364 


18.170 


13 


May 6 


592.8934 


1.19 


212 


17.841 


18.039 


18.182 


18.369 


18.165 


14 


May 11 


597.8990 


1.27 


192 


17.855 


18.013 


18.207 


18.369 


18.163 


15 


May 12 


598.9280 


0.84 


214 


17.836 


18.037 


18.175 


18.375 


18.160 


16 


May 13 


599.9053 


1.04 


167 


17.829 


18.022 


18.210 


18.374 


18.169 


17 


May 20 


606.9064 


0.90 


163 


17.848 


18.031 


18.208 


18.345 


18.156 


18 


May 21 


607.9159 


1.09 


182 


17.864 


18.014 


18.214 


18.373 


18.173 


19 


May 23 


609.9184 


1.33 


205 


17.857 


18.042 


18.205 


18.327 


18.161 


20 


May 26 


612.9050 


1.15 


220 


17.865 


18.059 


18.202 


18.346 


18.165 


21 


May 27 


613.9510 


1.15 


231 


17.853 


18.044 


18.181 


18.381 


18.164 


22 


May 28 


614.9155 


0.87 


238 


17.859 


18.059 


18.206 


18.352 


18.169 


23 


May 29 


616.0342 


1.35 


161 


17.848 


18.050 


18.215 


18.379 


18.171 


24 


Jun 5 


622.9486 


1.19 


238 


17.875 


18.073 


18.212 


18.359 


18.172 


25 


Jun 9 


626.9143 


1.36 


206 


17.918 


18.031 


18.229 


18.344 


18.168 


26 


Jun 11 


628.9096 


0.84 


196 


17.878 


18.041 


18.239 


18.363 


18.175 


27 


Jun 19 


636.9364 


0.90 


188 


17.887 


18.094 


18.266 


18.368 


18.176 


28 


Jun 23 


640.9468 


1.24 


228 


17.885 


18.069 


18.248 


18.380 


18.174 


29 


Jun 25 


642.9421 


1.29 


207 


17.922 


18.072 


18.256 


18.335 


18.160 


30 


Jul 5 


652.9233 


1.17 


232 


17.932 


18.076 


18.273 


18.401 


18.176 


31 


Jul 13 


660.9244 


1.06 


157 


17.950 


18.098 


18.302 


18.398 


18.178 


32 


Jul 21 


668.9129 


1.29 


164 


17.971 


18.067 


18.296 


18.403 


18.157 


33 


Jul 27 


674.9026 


1.27 


202 


17.915 


18.120 


18.288 


18.451 


18.176 



a All frames were taken in 2008. 

b FWHM of the seeing disc in arcsec. 

c SNR of the S40 field star within a circle of radius FWHM. 
d r— SDSS brightness of A in mag. The typical error is 0.010 mag. 
° r-SDSS brightness of B in mag. The typical error is 0.012 mag. 
f r-SDSS brightness of C in mag. The typical error is 0.018 mag. 
s r-SDSS brightness of D in mag. The typical error is 0.018 mag. 
h r-SDSS brightness of S40 in mag. The typical error is 0.006 mag. 
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TABLE 2 

Time delays (in days) of H1413+117. 



Method 


Simulations 








D 2 


NORMAL 


-17±| 


-19 ± 6 


24+5 




BOOTSTRAP 


-18 ± 3 


-23lf 


20 ± 5 


x 2 


NORMAL 


-17 ± 3 


-20 ± 4 


23 ± 4 




BOOTSTRAP 


-22+^ 


-23+.* 


21+ 4 



Note. — Ar;j = r,- - n, so B and C are leading, and D 
is trailing. All measurements are 68% confidence intervals. 
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TABLE 3 
Modelling results. 



Parameter MacLeod et al. (2009) This paper 



X 2 /dof 


4.9/5 


7.5/7 
0C68 


&G1 


C'66 


Aa G i 


-0'.'166 


-0'.'165 


AS G1 


0'.'556 


0'.'552 


eci 


0.26 


0.28 


Oeai 


-36? 5 


-37? 6 


&G2 


0'.'63 


0'.'45 




EE -l'/87 


EE -1"87 


A<5 G2 


EE 4'.' 14 


EE 4" 14 


le.xt 


0.087 


0.11 




50? 1 


45? 4 






1.88 



Note. — Position angles (ellipticity of Gl and 
external shear) are measured east of north and po- 
sitions are relative to image A (negative Aa values 
are eastward of image A). Here, b, e, 7, and zi de- 
note mass scale, ellipticity, shear strength, and lens 
rcdshift, respectively. 



